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Abstract 

We report on exclusive decays of B mesons into final states containing charmonium using data 
collected with the BABAR detector at the PEP-II storage rings. The charmonium states considered 
here are J/i/j , ■0(25'), and Xci- Branching fractions for several exclusive final states, a measurement 
of the decay amplitudes for the — > J/tp K* decay, and measurements of the B^ and B~^ masses 
are presented. All of the results we present here are preliminary. 
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1 Introduction 



An understanding of the decays of B mesons to final states including a charmonium resonance 
(J/-0, ijj{2S), Xci) is a prerequisite to an analysis of CP violation in the B system. In this paper 
we report the measurement of several branching fractions of exclusive decays, some of which have 
been used in our measurement of sin 2/3 ||T|. The channels considered are listed in Table |^. Here 
and throughout this paper the inclusion of charge conjugate states is implied. 



Table 1: B meson decay modes considered in this paper. 



Channel 



Secondary decay mode(s) 



B"^J/i;K^g 

B+^J/ipK+ 

B^^J/i^K*^ 

B+-^J/ijjK*+ 

B^^^{2S)Kl 

B+^il){2S)K+ 

B+^XciK+ 





e+r ; 












e+r ; 






£+£- ; 























We have used some of these exclusive modes to measure the masses of the B~^ and B^ mesons 
and their mass difference. We also present initial results on the yield of B^ J/ip which will 
be used for a future CP analysis. Finally we describe an amplitude analysis of the B J/ip K* 
decay. 

2 The BaBAR detector and dataset 

The BABAR detector is located at the PEP-II storage ring, an e~^e~ facility operating at the Stanford 
Linear Accelerator Center. PEP-II collides 9.0 GeV electrons with 3.1 GeV positrons to give a center 
of mass energy of 10.58 GeV, the mass of the T(45) resonance. 

The BABAR detector is described elsewhere [Q; here we give only a brief overview. Surrounding 
the interaction point is a 5-layer double-sided silicon vertex tracker (SVT) which gives precision 
spatial information for all charged particles, and is the primary detection device for low momentum 
charged particles. Outside the SVT, a 40-layer drift chamber (DCH) provides measurements of 
charged particle momenta. The dE/dx information from the DCH and SVT is used for particle 
identification. Beyond the outer radius of the DCH is a detector of internally reflected Cherenkov 
radiation (DIRC) which is used primarily for charged hadron identification. The detector consists 
of quartz bars in which Cherenkov light is produced as relativistic charged particles traverse the 
material. The light is internally reflected, and the Cherenkov rings are measured with an array 
of photo-multiplier tubes mounted on the rear of the detector. A CsI(Tl) crystal electro-magnetic 
calorimeter (EMC) is used to detect photons and neutral hadrons, as well as to identify electrons. 
The EMC is surrounded by a super-conducting solenoid which produces a 1.5 T magnetic field. The 
Instrumented Flux Return (IFR) consists of multiple layers of resistive plate chambers interleaved 
with the flux return iron. It is used in the identification of muons and neutral hadrons. 

The data used in these analyses correspond to a integrated luminosity of 7.7 fb~^ taken on the 
r(45) and 1.2 fb"^ taken 0.04 GeV below the peak. The data set contains 8.8 x 10*^ BB events. 
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For the analysis of the B meson masses we use a restricted set of data corresponding to 4.6 fb 

3 Particle reconstruction 

Inclusive charmonium reconstruction is described in detail in another contribution to this confer- 
ence Q. 

We here reconstruct J/ip candidates by combining pairs of oppositely charged tracks within the 
angular range 0.41 < 9 < 2.41(2.53) for electron (muon) candidates, where 9 is the polar angle to the 
beam axis. The invariant mass of the candidate must lie in the range 2.95 < mj/^ < 3.14GeV/c^ 
and 3.06 < mj/^ < 3.14GeV/c^ for decays to e+e" and /U"*"/!", respectively. When the J/tp 
decays to electrons, we demand that at least one of the tracks pass stringent particle identification 
requirements based on the ratio of the energy deposited in the EMC to the track momentum {E/p), 
and on the ionization loss of the track in the drift chamber (dE'/dx). For i-i~^fj.~ candidates, one 
track is required to pass a loose muon selection on the basis of the number of hit layers in the IFR 
and the other track is required to have an associated energy in the EMC which is consistent with 
a minimum ionizing particle. In the case of the decay to electrons, we apply a procedure to add 
photons which are close to the electron tracks and thereby reduce the impact of bremsstrahlung on 
the reconstruction efficiency. 

We select ^(25') candidal es in a similar way. For the decay to fi^fj, the invariant 

mass of the candidate is required to be within 0.05GeV/c^ of the nominal mass. In the case of 
■0(25) decays to e~^e~ the lower limit is relaxed to 0.25GeV/c^ below the nominal mass value. 
For the decay ip{2S) to J/jpTr~^TT~ , J/ip candidates are combined with pairs of oppositely charged 
tracks which originate from a common vertex, and the mass difference between the resulting ip{2S) 
candidate and the J/il) is required to be within 0.05 GeV/c^ of the nominal mass difference. For the 
decay J/i^^ e+e", the mass difference is relaxed to —0.25 < m^(2S') ~ < 0.05GeV/c^. 

Candidate Xci mesons are reconstructed via their decay to J/V'T- The 7 candidates are selected 
by requiring a neutral cluster in the EMC that has a distribution of crystal energies consistent with 
a 7 shower. The mass difference between the reconstructed Xci and the J/ij: is required to satisfy 
0.35 < AM < 0.45GeV/c^, and the momentum of the Xci in the T{AS) rest frame must lie in the 
range 1.15 < p* < 1.70GeV/c. 

Kg vr+vr^ candidates are formed from pairs of oppositely charged tracks which have an 
invariant mass between 0.489 and 0.507 GeV/c^. K'^ tt^tt^ candidates are required to have a 
mass between 0.470 and 0.525 GeV/c^ and an energy greater than 0.8 GeV. A vr*^ decay to two 
photons is observed in the EMC either as a single neutral cluster with substructure or as two 
distinct 7 clusters. The most probable decay point of the Kg is determined after refitting the two 
TT^ mesons at several points along the path defined by their summed momentum vector and the 
J/ip vertex. 

We reconstruct K*^' decays to K~^7r~ and Kgir^, and K*~^ decays to KgTT~^ and K~^tt^. In all 
cases the candidate K* is required to have an invariant mass within 0.075 GeV/ of the nominal 
value. 

A K^ candidate is reconstructed using neutral clusters observed in the EMC or the IFR. For 
an EMC candidate we require the deposited energy to be between 0.2 and 2.0 GeV and the cluster 
center-of-gravity to be well contained within the fiducial volume of the calorimeter (cos^ < 0.935). 
We reject candidates that are likely to be produced by photons by means of energy-dependent 
criteria based on the spatial distribution of the deposited energy in the cluster. A neutral cluster 
that can be combined with other neutral clusters to form a vr'^ or clusters with sub-structure 
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consistent with a tt^ is also rejected. A candidate observed in the IFR is required to have the 
cluster center within the fiducial volume (—0.75 < cos^ < 0.93) and to have a signal in at least 
two detector layers. We also apply additional isolation criteria to remove candidates that may have 
been split from clusters produced by charged particles. 

4 Exclusive B reconstruction 

For the decays 5°^J/^if°(7r+7r- ") and B'^—>-J/il^K'^, we require | cos^i^|, the absolute value of the 
cosine of the helicity angle of the J/ip , to be less than 0.8 (0.9) for J/ip decays to e+e" {p,~^n~). In 
addition we require that the Kg be consistent with having originated from the J/^p vertex. 

For the decay B^~^J/ipKg{TT^7r'^) we require \cos9h\ < 0.75 (0.8) for J/il^ decays to e+e~ 
{lJL^IJL~). To further reduce background in the J/if) ^ e^e~ channel, we require both tracks to 
satisfy electron identification criteria, one stringent and one loose. 

In reconstructing the decays B^—>-J/ipK*^ and B'^—fJ/tpK*'^ we require that a candidate 
charged kaon satisfy particle identification criteria based on ionization loss in the DCH and SVT 
and on the Cherenkov angle measured in the DIRC. The candidate 7r° mass is required to lie in the 
range 0.115 to 0.150 GeV/c^. We require | cos6t\ < 0.9, where 9t is the angle between the thrust 
direction of the reconstructed B and that of the rest of the event in the T(4S') rest frame. For the 
decay K*^ K^t:^ , the vertices of the K*^ and J/tp must be consistent with a single production 
point. 

For the decays B+^'i/j{2S)K+ and B'^^'ip{2S)K^{7r+n-) we require | cos 6't| < 0.9 for V(25) 
J/t/j7r~^7r~ decays and | cos 6h\ < 0.8 for ip{2S) decays. We also require the flight length 

to be greater than 2.5 mm and the to satisfy loose kaon identification criteria. 

We reconstruct the decay B'^^XciK'^ with the requirement |cos0t| < 0.9 and demand that 
the pass loose kaon identification criterion to further reduce the background in this channel. 

For B° JH^Kl decays we require | cos6h\ and | cos^^l to be less than 0.9, where 6b is the 
angle of the B candidate direction with respect to the beam axis in the rest-frame of the T(45'). 
We also require the sum of these quantities to be less than 1.3. 

To isolate the signal for each mode we use the variables Ai?, the difference between the re- 
constructed and expected B meson energy measured in the center-of-mass frame, and mES) the 
beam-energy substituted mass. These variables are defined as: 

mES = -Pf, (1) 

/^E = E%-El (2) 

where E^ is the beam energy in the center-of-mass, i.e., half the center-of-mass energy, and E"^ 
and p*B are the energy and momentum of the reconstructed B meson in the center-of-mass. For 
the B^ J/ijjK^ selection the momentum of the candidate is obtained by constraining the 
invariant mass of the and J/if) combination to the mass of the B meson. Therefore, in this mode 
only can be used to separate the signal from background. We exclude any event that passes 
the other exclusive B selections, has mES greater than 5.2GeV/c^ and |AE'[ less than 0.1 GeV. 

Only one exclusive candidate per event is accepted. If there are multiple candidates, we select 
the one with the smallest value of |A£^|. Exceptionally, in the B^ J/tpK^ selection we choose 
the candidate with the largest energy as measured by the EMC. If none of the candidate 
mesons have EMC information, we choose the candidate that has the largest number of layers with 
hits in the IFR. 
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5 Results 



5.1 Branching fraction measurements 

When deriving branching fractions we have used the secondary branching fractions and their asso- 
ciated errors pubhshed by the Particle Data Group 

We determine the number of BB events from the difference in the multi-hadron rate on and off 
the T{4S) resonance, normahzed to the respective luminosity. This leads to a systematic error of 
3.6% on all measured branching fractions. We have assumed the branching fraction of the T{4S) 
to BB is 100%, with an equal admixture of charged and neutral B final states. 

The efficiencies for each mode have been obtained from Monte Carlo simulations complemented 
with measurements of tracking and particle identification efficiencies extracted from the data. From 
particle identification control samples we assign a systematic error of 2% (3%) per electron (muon). 
We attribute a 5% systematic error to the vr*^ reconstruction efficiency and resolution. The track 
finding efficiency has an uncertainty of 2.5% per track. Uncertainties in the modeling of the track 
resolution lead to an additional 1-2% error depending on the details of the primary and secondary 
decays. 

For each mode the shape of the beam-energy substituted mass distribution is parameterized 
with the sum of a Gaussian and the ARGUS function |^] . We assign a systematic error due to our 
uncertainty of the shape of the background of between 1% and 9% depending on the mode. For the 
B — > J/tpK* channels, a likelihood fit is performed for all the decay modes simultaneously, taking 
into account the cross-feed between decays. 

Figure |l| shows the ttt-es and AE distributions of the candidates. In Table ^ we present the 
yields and measured branching fractions for the individual exclusive modes. Figure § shows the 
measured branching fractions compared to the values compiled by the Particle Data Group 0]. 

Table 2: The yields and measured branching fractions for exclusive decays of B mesons involving 
charmonium. The yield only includes the statistical error. For the branching fractions, the first 
error is statistical and the second systematic. All results are preliminary. 



Channel 




Yield 


Branching fraction/10 ^ 


B^-^J/tpK" 




^ vr+TT- 93 ± 10 


10.2 ± 1.1 ± 1.3 






^ 7r°7r0 14 ± 4 


7.5 ± 2.0 ± 1.2 


B+^J/'4)K+ 




445 ± 21 


11.2 ± 0.5 ± 1.1 


B^^J/ijK*^ 




188 ± 14 


13.8 ±1.1 ±1.8 


B+^J/4^K*+ 




126 ± 12 


13.2 ±1.4 ±2.1 


B°^tP{2S)K^ 




23 ± 5 


8.8 ± 1.9 ± 1.8 


B+^'il^{2S)K+ 




73 ± 8 


6.3 ± 0.7 ± 1.2 


B+^XciK+ 




44 ± 9 


7.7 ± 1.6 ±0.9 



5.2 Observation of a signal for the decay — > J/%I^K^ 

Figure |3| shows the /S.E distributions for the B^ J/ipK'^ candidates in data and Monte Carlo 
events. We determine the yield by counting events with AE < 0.01 GeV and subtracting the 
background contributions. There are two categories of background to this mode. The first arises 
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Figure 1: Distributions of candidate events in rriEs and AS. The upper plot shows the modes 
and the lower plot the modes. All results are preliminary. 
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Figure 2: Summary of branching fraction measurements for charmonium + K channels and 
comparsons with the PDG 2000 values. All results are preliminary. 
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Figure 3: Comparison of data with Monte Carlo simulation for the —>■ J/ijjK^ selection. The 
results presented here are preliminary. 



from other B J/'ipX decays and is estimated from simulations normalized to the measured 
J/'0 yield. We estimate the systematic error on the yield from this background determination by 
varying the reconstruction efficiency, the branching fractions of the major background modes, 
and the helicity amplitudes used in the simulation of i? — > J/i/jK* decays. The second background 
arises from non-J/ip modes and is measured from a sideband above the J/^lJ peak in the data. The 
systematic error from this contribution is determined by varying the shape of the background. We 
measure the yield of B^ — > J/tpK^ events to be 82 it 14 (stat) ±9 (syst). This is in good agreement 
with the expected number of 93 signal events predicted from Monte Carlo simulation. 

5.3 Measurements of the B° and masses 

Measurements of the B^ and B^ invariant masses have been performed using the decay modes 
S°^J/'0K°(7r+7r-), B°^J/ijK*^{K+7r-) and B+^J/ijK+. These modes are chosen because 
they have small backgrounds and the masses of the secondary decay products are well known. The 
event selections are as described in section 5.1, with the additional requirements niEs > 5.27GeV/c^ 
and \AE\ < 36 MeV. 

The B candidate invariant mass is derived by fitting the decay products to a common vertex, 
with the masses of the J/^p and Kg constrained to their nominal values. Uncertainties in the 
magnetic field and in the internal and relative alignment of the tracking devices can introduce a bias 
in the momentum measurement. The size of this effect is quantified by comparing the reconstructed 
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mass of J/ip^ A*^^^ candidates, determined by fitting to the invariant mass distribution, to the 
nominal mass value. Any observed shift is subsequently applied to the track momenta in simulated 
data to determine a correction to the measured B mass. The systematic error attributed to this 
correction is derived from statistical uncertainty on the parameters from the fit to the J/ij) invariant 
mass. The i?*^— > J/-0i^g(7r"'"7r~) sample requires special consideration as the decay products of the 

do not come from the interaction point and are sensitive to details of the track parameterization. 
A fit to the observed tt+vt^ invariant mass is performed and a correction derived in the same way 
as described above. The correction applied to the B^—fJ/ip Kg{TT^7r~) sample is taken to be the 
mean of the correction factors determined from the fits to the J/ip and distributions, taking the 
semi-dispersion as the error. The resulting systematic uncertainty is ±0.62, 11q^ and j]g44MeV/c^ 
for the B'^^J/i;K^{TT+TT-), B^^J/4^ K*'^{K+-k-) and B+^J/'ipK+ modes, respectively.' 

An additional uncertainty comes from background contamination in the event samples, which 
are determined from a fit to the sideband events, and is found to be between 2% and 4%. The 
measurement of the mass is performed by fitting a single Gaussian and a flat background to the 
B invariant mass distribution. Examples of the mass distributions are shown in Fig. ^. The 
distortion of the mass measurement due to the presence of background has been estimated by 
removing the events with smallest mass and the N events with highest mass, where N is the 
number of background events in the sample. The corresponding systematic uncertainty is ~^q'qq, 
toil and toil MeV/c2 for the J/^ K0(7r+^- 
respectively. 



B'^^J/ijK*^{K+TT') and B+^J/il)K+ modes. 
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Figure 4: The reconstructed B mass distribution for the (a) B~^^J/ip , (b) 

B^^J/ip K*^{K^-K~) and (c) B^^J/ip K'^{-k^-k~) samples. The results presented here are pre- 
liminary. 



The B masses have been measured to be: 

m(5°) = 5279.0 ± 0.8 MeV/c^ 
m{B+) = 5278.8 ± 0.6 1[5:^ MeV/c^, 
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where the first error is the quadratic sum of the statistical and uncorrelated systematic errors and 
the second error is the correlated systematic error. 

The mass difference between the and mesons is evaluated by fitting the tues distributions 
of the three above-mentioned channels with the ARGUS function to describe the background and 
a Gaussian to describe the signal. The use of m-ES has the advantage that it reduces the sensitivity 
to the measured momentum scale and the uncertainties in the energy scale of the beam particles 
cancel in the mass difference. The same systematic study has been performed as in the invariant 
mass measurement described above, and was found to contribute 0.01 MeV/c^ to the systematic 
error on the mass difference measurement. Simulations indicate that the effect of the uncertainties 
in the beam parameters on the mass difference measurement is only 0.001 MeV/c^. 

We also consider how the uncertainty in the shape of the background under the signal affects the 
mass difference measurement. We estimate the uncertainty by fitting the shape of the distribution 
in the AE sidebands and using these parameters when fitting to the signal. The effect on the mass 
difference between fixing the background shape or not is found to be 0.04MeV/c^. 

We measure the mass difference to be: 

m(5°) - m(B+) = 0.28 ± 0.21 ± 0.04 MeV/c^ 

where the first error is the quadratic sum of the statistical and uncorrelated systematic errors and 
the second error is the correlated systematic error. 

5.4 Angular analysis of S — ^ J/iJ^K* 

The B —5- J/ipK* decay proceeds through three amplitudes, corresponding to the three different 
helicity configurations of the decay products . The transversity formalism involves linear combi- 
nations of these amplitudes, denoted by Aq, At and ^y. Both and are CP even while At is 
CP odd. The size of the CP odd contribution in the decay must be known before a value of sin 2/3 
can be extracted from this decay channel. 

The event selection is similar to that used for the branching fraction measurement in section |5.1| . 
In this analysis we have considered only those channels which have a final state composed solely of 
charged particles. B candidates were required to have a reconstructed mass within ibO.OlGeV/c^ 
of the nominal value and \AE\ < 0.075 GeV. 

The background is determined from the sideband region m^s < 5.25 GeV/ c^. The amplitudes 
are determined from a fit using the unbinned extended likelihood method 0, that takes into 
account the normalization condition |^oP + + l^tP = 1; the finite detector acceptance and 
the background contributions (assumed to have a flat angular distribution). 

The 68% contours of the fit are presented in the |AoP + |^||P + |^tP = l plane of the |^oPj 
l^ltp space in Fig. ^ The fraction of the amplitude that is CP odd is determined to be 
i^tp = 0.13 ± 0.06 ± 0.02, while the longitudinal polarization (Tl/T) is found to be |^oP = 
0.60 ±0.06 ±0.04. 

Systematic errors arising from our knowledge of the background, the acceptance corrections, 
the cross- feed among B J/tpK* modes, and the contribution due to heavier K* mesons have 
been considered. The effect on the transversity amplitudes and phases are summarized in Table |3|. 

The results obtained are presented in Table ^. 
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Figure 5: Transversity angle fits. The modulii of the amplitudes are shown as 68% contours in the 
top plot while the relative phase of the amplitudes are displayed in the bottom plot. Also shown 
are the results from CLEO [Hi. 



Table 3: Systematic uncertainties in the measurement of transversity amplitudes. 



Source 


|ylo|^ l^lll^ V^ll (rad) (rad) 


Monte Carlo statistics 
Backgrounds 
Angular acceptance 
Cross-feed background 
Heavy K* 


0.014 0.014 0.016 0.12 0.08 
0.011 0.009 0.001 0.01 0.03 
0.020 0.011 0.020 0.13 0.05 
0.025 0.006 0.030 0.02 0.03 
0.011 0.004 0.007 0.07 0.01 


Total 


0.038 0.021 0.040 0.19 0.10 



Table 4: Measured transversity amplitudes. The third amplitude is determined from the normal- 
ization condition while the phase 0o is, by convention, set to zero. 







(rad) 


ipt (rad) 


0.60 ± 0.06 ± 0.04 


0.13 ±0.06 ±0.02 


2.58 ±0.39 ±0.20 


0.01 ±0.27 ±0.10 



6 Summary 

We have presented preliminary measurements of branching fractions of B mesons to several two 
body charmonium final states. The results are in good agreement with previous measurements. A 
signal for the decay Jj'\\)K^ is observed, with a yield compatible with our expectation. 

The 5° and S+ masses and their mass difference have been measured, with results that are in 
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good agreement with the world average values 

Finally we have presented an analysis of the transversity amplitudes in the decay B J/tp K* 
that confirm that this final state is dominantly CP even and the CP asymmetry measurements in 
this channel will have a small dilution. 
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